Anomalous baryon production in heavy ion collisions at LHC energies 
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In ultrarelativistic heavy ion collisions the baryon-meson ratio is not fully understood in the 
intermediate-to-high transverse momentum region. Although quark coalescence methods combined 
with jet fragmentation yield results close to the experimental data, some tendencies are not repro- 
duced and require further investigation. We propose a new channel, namely extra quark-antiquark 
and diquark-antidiquark pair production from coherent gluon field created in the early stage of heavy 
ion collisions. This process is described quantitatively by the time-dependent Schwinger-mechanism. 
The extra diquark-quark coalescence yields are able to explain the measured anomalous baryon pro- 
duction. 



I. INTRODUCTION 

The aim of ongoing investigations at ultrarelativistic 
heavy ion colliders e.g. the BNL Relativistic Heavy Ion 
Collider (RHIC), and the CERN Large Hadron Collider 
(LHC), with up to y/s = 200 AGeV and = 5.5 ATeV 
center of mass energies respectively, is to create an ex- 
tremely high energy density, comparable to that of the 
early Universe. At such high energies the colliding nu- 
clei can be modeled by two strongly Lorentz-contracted 
sheets of nucleons, with 7 « 100 at RHIC and 7 2750 
at LHC, surrounded by a gluon cloud. Their overlap cre- 
ates a strong and rapidly changing color field. Particle 
(eg. quark-antiquark) pairs are produced in this field, a 
process similar to the Schwinger-mechanism in Quantum 
Electrodynamics (QED). 

The phenomenology of particle production in high en- 
ergy pp and heavy ion collisions traditionally relies on 
the introduction of Quantum Chromodynamics (QCD) 
strings or ropes. The string picture provides a mech- 
anism to convert the collisional energy to field energy, 
which in turn is converted to particle pairs. Simple mod- 
els 0, Q are based on a time-independent picture and 
work properly at ISR and SPS energies. However at ex- 
treme high energies, where the field amplitude is rapidly 
changing, like at LHC, this is not a realistic assumption. 

After creation, the particles start to interact with the 
surrounding matter, so low-momenta light and strange 
quarks lose most of their original momentum, and ther- 
malize. In contrast, high-momentum particles quickly 
escape from the strongly interacting matter, and main- 
tain their original momentum spectra. This way, the 
measured high momentum hadron spectra will be closely 
related to that of the primordial quarks. 

However, in a strong field, not only quark-antiquark 
pairs can be created, but composite objects, such as di- 
quarks too. We are focusing here on primordial diquark- 
antidiquark yields, as an extra source to final state 
hadron production at intermediate transverse momenta. 

We build a simple phenomenological framework for the 
description of mid-to-high transverse momenta hadron 
production in heavy ion collisions, focusing on LHC 



energies. In addition to the usual thermal coales- 
cence/recombination|4j49| and hard scattering (pertur- 
bative QCD) @, OJ, M,[M , 

we consider new coalescence 
channels, that are based on a dynamical model of quark 
and diquark pair production in strong fields [37l l4l[ |42| . 

After reproducing the existing experimental data from 
the ALICE and CMS experiments on unidentified hadron 
yields, we argue, that the inclusion of these new chan- 
nels can lead to anomalous (enhanced) baryon-to-meson 
ratios in the mid-to-high transverse momentum range. 
Such anomalous ratios have already been observed at 
RHIC [l7l - [2l| . From our model, we expect a similar ef- 
fect at LHC, where the experimental data on identified 
hadron spectra at such momenta is not yet available, but 
will be, after the accumulation of sufficient statistics. 

This paper is structured in the following way. In 
the 2nd section we briefly review the usual description 
of hadron production, including the thermal coalescence 
and the jet fragmentation channels, and show our results 
obtained with them. The 3rd section reviews the tra- 
ditional string phenomenology, followed by the Wigner- 
function based derivation of the time-dependent model 
of pair production. In the 4th section we present the 
numerical results obtained by our framework. We close 
the paper with a brief discussion of our results, and an 
outlook to possible future investigations. 



II. PHENOMENOLOGY OF HADRON 
PRODUCTION AT MID-TO-HIGH p T 

In proton-proton collisions, the hadron produc- 
tion is usually described by string fragmentation at 
Pr < 2 GeV/c and jet fragmentation at pr > 2 GeV/c 
[3j. In heavy ion collisions a large number of the produced 
partons are thermalized, creating the wanted quark-gluon 
plasma. These low momentum thermal quarks can form 
hadrons with higher momentum. This process is usually 
called quark coalescence 0d3- This channel will dom- 
inate the medium transverse momentum spectra. The 
validity of this picture was supported by RHIC data on 
proton-to-pion ratios [l7l - [2]] | and the elliptic flow of iden- 
tified hadrons [22[ ■ Another consequence of this dense 
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partem matter is the appearance of jet energy loss [Hj]. 

In this section, we briefly review these aspects of 
hadron production in heavy ion collisions, especially at 
LHC energies. 

A. Jet fragmentation 

In proton-proton collisions, the production of the most 
abundant hadrons can be calculated by the QCD im- 
proved parton model 0, El, 13 • At first, we determine 
the primary quark and gluon distributions. Later on, 



We use LO MSTW parton distribution function [l5| 
and AKK fragmentation functions [l6[ . In heavy ion col- 
lisions Ta(6) — J dzpA(b, z) is the nuclear thickness func- 
tion normalized as usual / d 2 6T/i(&) = A. The Woods- 
Saxon formula is applied for the nuclear density function 
PA(b, z). The integral in b indicates the nuclear overlap 
and the consideration of the Glauber geometry. Focus- 
ing on high momentum parton production, we neglect 
nuclear shadowing, so the parton distribution of the nu- 
clei reads: 

fa/A{x) = (^fa/p(x) + ^fa/n{x)j ■ (3) 

The average jet energy loss is taken into account with 
the simplified GLV formula (l2| : 

where Cp is the Casimir-factor of the jet, N(E) is a 
smooth function of energy, calculable with the help of 
(l2^ . A is the radiated gluon free path, jj 2 /A oc a 2 p is a 
transport coefficient of the medium proportional to the 
parton density p, and L is the plasma thickness. We 
concentrate on midrapidity {y = 0), where the jet trans- 
verse momentum is shifted before fragmentation by the 
medium to p* = pc — A_E, which changes the z c param- 
eter in the integrand to z* = z c j [\ — ■ Then the 

hard reaction scale is Q = £p-r and the fragmentation 
scale is Qp = Q/z*. The parameter £ will be fitted to 
match the pp data, and the same value is used for the 
heavy ion calculations. 

By means of this model, we can describe the very high 



we fragment these high energy partons by means of frag- 
mentation functions. In heavy ion collisions an additional 
process, namely jet energy loss needs to be included due 
to the produced hot dense deconfined matter. 

In our calculation, we use a leading order pQCD frame- 
work, to make the transition to heavy ion collisions sim- 
pler. In A+B heavy ion collisions, we consider a geo- 
metrical superposition of nucleon-nucleon collisions, and 
include the effect of jet energy loss by shifting the frag- 
mentation function. 

The basic equations for calculating hadron spectra 
read: 



fi), (!) 



(2) 



I 

momentum hadron production from jets. Our numerical 
results will be described later. 



B. Quark coalescence 

In the heavy ion collisions a dense parton matter 
will be produced. This matter will thermalize and ex- 
pand, reaching the deconfinement transition tempera- 
ture known to be T w 180 MeV, calculated on the 
lattice [l(| EH- At this point, the combination of par- 
ticipant quarks results in the dominant hadron produc- 
tion yield at low momenta, and give a significant contri- 
bution in the intermediate momentum region, overlap- 
ping with jet fragmentation. This yield will be deter- 
mined by quark coalescence 0-0 applied in the low-to- 
intermediate transverse momentum range. Among the 
numerous models, we will use a simple, phase-space dis- 
tribution based description. We assume an infinitely 
sharp (pre-)hadron momentum distribution, so the cal- 
culated meson and baryon spectra reads: 

s W = (f^X P ^ /a(r '^ /2)/b(r ^ /2) ' 

(5) 

E W = jtri P ^/a(r, P/3)/ 6 (r, P/3)/ c (r, P/3), 

(6) 

where gm and gp are the degeneracies of the resulting 
hadron, g n + = 1, g p + = 2, g p + = 3, gj+ = 3 etc. The 
spacelike hypersurface S over which we integrate is cho- 
sen to be that of constant proper time r = \Jt 2 — z 2 , and 



Pfe (pp-Hi + X) = EabcdI^adx b dzJ a/p (x a ,Q 2 )f b/p (x b ,Q 2 )D h/c (z c ,Q 2 F )^ ~* C S(s + t + 

E^(AB^h + X) = 

= J d 2 b d 2 rT A (r)T B (\b - r\) £ Qbcd / dx a dx b f a/A (x a , Q 2 )f b/B (x bl Q^D h/c (z*,Q%)^ (^Y^ 
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FIG. 1. (Color online.) Measured charged hadron transverse 
momentum spectra at 2.76 ATeV in central PbPb collisions 
[2cj |. Calculation includes quark coalescence and parton frag- 
mentation yields with jet energy loss from GLV model (l2l | 
with opacity L/X = 5. 
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FIG. 2. (Color online.) Nuclear modification factor (Raa) in 
a logarithmic scale for charged hadrons, pions and protons ac- 
cording to our latest calculations including quark coalescence 
and jet energy loss in central PbPb collisions at 2.76 ATeV. 
Data are from ALICE HHHl. 
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FIG. 3. (Color online.) Nuclear modification factor (Raa) in 
a linear scale, obtained from simple perturbative-QCD calcu- 
lations including jet energy loss in central PbPb collisions at 
2.76 ATeV, where the experimental data are from ALICE [33j 
and CMS \M- 



the one parton Wigner-functions are assumed to have the 
exponential form: 



fi(r,p) = ll e-v v / T 0( Pa -p)e-^ 



(7) 



where the flow four-velocity v follows a Bjorken sce- 
nario: 



( coshf? cosh^T, 
sinh r\x cos <fi, 
amhrjT sin</>, 
sinh?y cosher ) 



(8) 



In this case the integral is completely insensitive to the 
parameter A. Essentially our free parameters will be the 



t/9q7t = tAt coalescence volume, and the vx transverse 
flow velocity given by tanhr^ = vx- These two can be 
fitted to the Raa of unidentified hadrons, giving a pre- 
diction for the identified spectra. 



C. The interplay of pQCD and coalescence 

Hadron production, especially charged pion, proton 
and antiproton production have been measured with 
high precision at RHIC energies in pp and AuAu col- 
lisions The analysis of the measured proton- 
to-pion ratio proved the unsatisfactory performance of 
perturbative QCD and independent jet fragmentation in 
the transverse momentum region of 2 GeV/c < px < 
8 GcV/c (23[. The introduction of parton coales- 
cence/recombination models [IHII offered a possible so- 
lution and these calculations displayed the dominance of 
this channel in the intermediate-py region. The result- 
ing yields of our model including fragmentation and co- 
alescence is shown on Figure 1. Figure 2. displays the 
nuclear modification factor, Raa- Figure 3 shows the ex- 
treme high momentum behavior of Raa obtained from a 
pure pQCD calculation. 

However, the data on Raa at RHIC energies (20l . |2~0 | 
displayed another anomaly: the Raa(p + ) > Raa{^ + ) in 
the window of 2 GeV/c < px < 8 GeV/c and even be- 
yond. This result contradicts the expectations from jet 
energy loss descriptions (see e.g. the GLV-model 112ft') 
and it was explained by quark-gluon jet conversion [24 . 
I2IH and hadrochemistry inside jet-cones J2(|. Further- 
more, a fine-tuned coalescence/recombination calculation 
may reproduce these data, however, we would expect that 
this channel disappears at px > 8 GeV/c (27|. With the 
appearance of new ALICE data at LHC energies (28| . 
the question of nuclear modification factors for identified 
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FIG. 4. {Color online.) The external field used for the pair 
production channel, see eq. [T8l 
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FIG. 5. [Color online.) Particle spectra from the Schwinger- 
mechanism for quarks and diquarks. 



hadrons appears again and it demands a detailed analy- 
sis. In the following section we propose another process 
that may contribute to the aforementioned anomaly. 



III. PAIR PRODUCTION IN STRONG FIELDS 

A. Constant fields and string phenomenology 

Early models of hadron production in proton-proton 
collisions, e.g. the LUND and FRITIOF models 0, Q, 
were based on q — q and qq — qq string formation. 
These strings are broken by the production of new 
quark-antiquark pairs. These models suppose a time- 
independent color field, resulting in the q — q and qq — qq 
yields: 



d 7 N 



dtd 3 xd 3 p 



2 exp I — 



Tr(m 2 +Pt) 



(9) 



In this model, characterized by the k string constant, 
both the diquark and high momentum quark production 
are suppressed. 

The description of the early stage of heavy ion col- 
lisions is different from the string formation in proton- 
proton collision, because of the higher density. For a 
very short time an extremely strong and rapidly changing 
color field will manifest. This modification was mimicked 
by introducing color ropes |35[ , still described by equa- 
tion ©, but characterized by a larger string constant 
(k ^> 1 GeV/fm). However this will not change the basic 
suppression behavior. Recently, it has been shown, that 
introducing a flavour dependent string (rope) values the 
measured baryon-meson yields can be approximately re- 
produced [43[. We choose a different approach, relaxing 
the assumption of a constant field. 



B. Time-dependent fields 

In this study, we use the Quantum Kinetic (QK) model 
for the description of pair production. It was shown, 



that certain time-dependent external fields can result in 
power-law spectra [37| . The QK method can be extended 
to non-Abelian cases, but it was shown, that the be- 
haviour has strong Abelian dominance (38|. This is the 
basis of modeling the quark and diquark production from 
strong SU(3) fields with QED based calculations, because 
the spectra is expected to be the same except for the am- 
plitude scale. 

The quantum kinetic equation for the non-Abelian case 
can be derived from the time evolution of the correspond- 
ing Wigner function W(k,t). The non-localities can be 
avoided by applying the gradient approximation (see (38j 
for details). The time evolution is then described by: 



d t W 



9 d 



(4{W, F 0i }- 



8dki 

+2{i^,[W )7 V]} - [F W ,{W/ 7 V}]) = 

= ifefrV, W} - im[ 7 °, W] + ig [A, , [ 7 V , W}] . (10) 

were m denotes the mass of the particles in considera- 
tion, g is the coupling constant, is the 4-potential 
of the external space-homogeneous color field and F^ is 
corresponding field tensor: 



d ll A v - dyA^ - iglAp, A v \. 



(11) 



While it is convenient to perform spinor and color de- 
composition of the Wigner function for direct numerical 
computations, if the symmetry of the external field is ex- 
ploited, the kinetic equations can be recast in the form of 
the U(l) case, where the final equations for the particle 
density / read [39j : 



dt 

dv 

dt 
da 

~dt 



e£e± 
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Here u and v are auxiliary functions and 



UJ 2 (q,t) 



(12) 
(13) 
(14) 

(15) 
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FIG. 6. (Color online.) Anomalous proton and pion yields 
from strong color fields, comparing to pQCD yields and AL- 
ICE data 



FIG. 7. (Color online.) The modification of the proton/pion 
ratio at LHC energies after introducing the anomalous proton 
(and pion) yield connected to pair creation from strong field. 



e\ = m 2 + ql 



(16) 



IV. NUMERICAL RESULTS 



q = ipx,p\\ - eA(t)) 



(17) 



The initial conditions at t — — oo are f(q) = v(q) = 
u(q) = 0. The only input for the QK model is the elec- 
tric held £(t) and its vectorpotential A(t). We chose the 
following time dependence (see Fig. 



£(t) 



^ 2 ± 



Ea[l 



t < 



t > 0. 



(18) 



where we use the phenomenological field parameters (E , 
(5, T2) according to our expectations: the first part of 
the formula is motivated by the increasingly overlapping 
nucleons, avoiding a full self-consistent QCD description. 
The decay of the field is suggested by Bjorken hydrody- 
namics. The time-scale of the raise and decay is n and T2 
respectively. The field amplitude is measured in critical 
Schwinger field units: E cr = m 2 /g, set by the effective 
quark mass. We chose m q = 300 GeV. The QK equa- 
tions are integrated numerically in time and in the longi- 
tudinal direction, giving f(pr) at t = +00 for quark and 
diquarks. Figure [3] illustrates the obtained q and qq spec- 
tra. Note, that at high momenta, the diquark and quark 
spectra are similar, and both are enhanced compared to 
the case of a constant field. 

The input of the coalescence equations ([5]) will then 
be: 



fs(r,p) = f(p T )Q{po - p) exp 



2A 2 



(19) 



where the parameters po and A are the same as in the 
calculation of the thermal spectra, with the addition of a 
diquark + quark — > hadron channel, governed by equa- 
tions formally equivalent to quark + quark — > meson co- 
alescence, only with a quark distribution interchanged 
with the diquark one. 



Earlier studies of particle production from time- 
dependent strong coherent fields [H, |4l| indicated that 
in case of a rapidly changing field, the obtained spectra 
of quark-antiquark pairs could easily become power law, 
instead of the Gaussian drop of usual Schwinger spec- 
tra from static strong electric field [43[. Thus we have 
an extra channel of quarks and antiquarks with momen- 
tum spectra different from usual pQCD results. Further- 
more, diquark-antidiquark pairs can be created, open- 
ing an anomalous channel for baryon production. Since 
quark and diquark yield are very close to each other in 
a time dependent field, the proton/pion ratio would be 
close to unity for hadron coalescence from pair-creation 
channels. Figure[5]displays these effects. This anomalous 
yield of protons from diquarks can overwhelm the frag- 
mentation component, as can be seen on Figs. [5] and [51 
where only the pQCD and pair creation channels are in- 
cluded. The parameters of the time dependence (n, T2, S) 
and the strength of the electric field (£0), wn l determine 
the resulting quark and diquark spectra. In our fit on 
the unidetified hadron Raa [!H we found Eq = 2.2E cr , 
d = 6.0, Ti = 0.31/m 9 , T2 = 0.36/m g for the external 
field, and tA t = 4226 fm 3 and v T = 0.7 for the thermal 
coalescence parameters. 

Finally, Fig. [5] displays Raa obtained from the com- 
bined framework, including parton fragmentation and 
quark coalescence using thermal source and pair produc- 
tion in time-dependent strong fields. In this way we can 
approximately reproduce the measured nuclear modifi- 
cation factor in central PbPb collisions at 2.76 ATcV, 
and a modified proton and pion suppression pattern will 
appear, namely Raa (proton) > RAA(pion), as shown on 
Fig. [5]) . Naturally, this anomalous proton yield will mod- 
ify the proton-to-pion ratio, although the modification 
could remain relatively small, as indicated by Figure [7]). 
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FIG. 8. (Color online.) The modification of the nuclear mod- 
ification factor (Raa) connected to the appearance of the 
anomalous hadron yields in PbPb collisions at 2.76 ATeV, 
compared to ALICE data [2§ |. 



FIG. 9. (Color online.) The nuclear modification factor 
(Raa) after including jet fragmentation, quark coalescence 
and the yield from time-dependent strong coherent field. The 
theoretical results are compared to ALICE data [28j|. 



V. CONCLUSIONS 



We have studied hadron production in heavy-ion col- 
lisions focusing on baryon-to-meson ratios and proton, 
pion Raa- We used a model combining jet fragmenta- 
tion with jet quenching and coalescence in a thermal 
quark matter characterized by the deconfinement tran- 
sition temperature T ss 180 MeV. We also introduced 
a new channel into the coalescence model, originating 
from extra quark and diquark yields from strong time- 
dependent fields. We have demonstrated that such a 
model can reproduce existing experimental data at LHC 
energies, and that an anomalous production, formerly 
seen at RHIC [2(| HH, could also be seen in future LHC 
data. 

High precision measurements of identified charged 



hadrons at LHC energies could become the basis of a 
proper discussion on the formation of coherent field and 
the existence of anomalous quark and diquark production 
channels. 

Future numerical studies could include a more realistic 
time and space dependence of the external gluon field 
(see e.g. (45l . |46| ) which may modify the spectra of 
produced quarks and diquarks. A similar investigation 
was discussed for the Abelian case in Ref. fiij . 
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